Commensurability effects in superconducting Nb films with quasiperiodic pinning 

arrays 
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We study experimentally the critical depinning current I c versus applied magnetic field B in Nb 
thin films which contain 2D arrays of circular antidots placed on the nodes of quasiperiodic (QP) 
fivefold Penrose lattices. Close to the transition temperature T c we observe matching of the vortex 
lattice with the QP pinning array, confirming essential features in the Ic(B) patterns as predicted 
by Misko et al. [Phys. Rev. Lett. 95(2005)]. We find a significant enhancement in I C (B) for QP 
pinning arrays in comparison to I c in samples with randomly distributed antidots or no antidots. 

PACS numbers: 74.25.Qt, 74.25.Sv, 74.70.Ad, 74.78.Na 



The formation of Abrikosov vortices in the mixed state 
of type-II superconductors Q and their arrangement in 
various types of "vortex-phases", ranging from the or- 
dered, triangular Abrikosov lattice to disordered phases 
S 3 has a strong impact on the electric properties of 
superconductors. Both, in terms of device applications 
and with respect to the fundamental physical properties 
of so-called "vortex-matter", the interaction of vortices 
with defects, which act as pinning sites, plays an impor- 
tant role. Recent progress in the fabrication of nanostruc- 
tures provided the possibility to realize superconducting 
thin films which contain artificial defects as pinning sites 
with well-defined size, geometry and spatial arrangement. 
In particular, artificially produced periodic arrays of sub- 
micron holes (antidots) |E IE S @ an d magnetic dots 
[EE! IHHil as pinning sites have been intensively inves- 
tigated during the last years, to address the fundamental 
question how vortex pinning - and thus the critical cur- 
rent density j c in superconductors - can be drastically 
increased. 

In this context, it has been shown that a very stable 
vortex configuration, and hence an enhancement of the 
critical current I c occurs when the vortex lattice is com- 
mensurate with the underlying periodic pinning array. 
This situation occurs in particular at the so-called first 
matching field Bi = $>o/A, i.e., when the applied field B 
corresponds to one flux quantum $0 = h/2e per unit-cell 
area A of the pinning array. In general, I C (B) may show 
a strongly non-monotonic behavior, with local maxima 
at matching fields B m = mB\ (m: integer or a rational 
number), which reflects the periodicity of the array of 
artificial pinning sites. 

As pointed out by Misko et al. ^| , an enhancement of 
I c occurs only for an applied field close to matching fields, 
which makes it desirable to use artificial pinning arrays 
with many built-in periods, in order to provide either 
very many peaks in I C {B) or an extremely broad peak in 
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I C (B). Accordingly, Misko et al. studied analytically and 
by numerical simulation vortex pinning by quasiperiodic 
chains and by 2D pinning arrays, the latter forming a 
fivefold Penrose lattice [lj], and they predicted that a 
Penrose lattice of pinning sites can provide an enormous 
enhancement of I c , even compared to triangular and ran- 
dom pinning arrays. 

We note that the discovery of quasicrystals has, 
until today, stimulated intensive investigations of a large 
variety of artificially generated quasiperiodic systems, 
such as semiconductor heterostructures |16I. optical su- 
perlattices 0] , photonic quasicrystals |l8j. atoms in op- 
tical potentials 0]i superconducting wire networks [2£j 
and Josephson junction arrays [2l| • The investigation of 
the static and dynamic properties of vortex quasicrystals, 
including phase transitions which may be tuned by tem- 
perature and magnetic field, is interesting, both from a 
practical point of view (regarding controllability and en- 
hancement of critical currents in superconductors) and 
also with respect to our understanding of fundamental 
aspects related to the physics of quasicrystals. 

Here, we present results on the experimental investi- 
gation of matching effects in superconducting Nb thin 
films containing various types of antidot configurations. 
We studied I C (B) at variable temperature T close to the 
superconducting transition temperature T c , and we com- 
pare Penrose lattices with triangular lattices, with ran- 
dom arrangements of antidots and with thin films with- 
out antidots. Our experimental results on Penrose ar- 
rays confirm essential features in the I C {B) patterns as 
predicted in . 

The experiments were carried out on d — 60 nm thick 
Nb films which were deposited by dc magnetron sput- 
tering in the same run on five separate substrates. Pat- 
terning was performed by e-beam lithography and lift-off 
to produce cross-shaped Nb bridges with circular anti- 
dots arranged in different geometries. Fig. QJa) shows 
the geometry of the Nb bridges, with segments of width 
w = 200 /im and length (separation between voltage 
pads) I = 600 ^im. On each chip we are able to directly 
compare eight different cross structures. Six of them con- 
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FIG. 1: (Color online) Images of (a) cross-shaped Nb bridge 
and of different arrangements of antidots (125 nm radius): tri- 
angular lattice (b); Penrose lattice (c); random arrangement 
(d). The lines in (b),(c) illustrate the lattice geometries. 



tain approximately N p =110,000 circular antidots with 
radius r — 125 nm or 200 nm, arranged in either a trian- 
gular lattice [Fig.^b)], a Penrose lattice [Fig.^c)], or in 
a random arrangement [Fig.^d)]. All samples have the 
same average antidot density n p = 0.52 /im" 2 , which cor- 
responds to a first matching field B\ = n p $o = 1.08 mT. 
For reference measurements each chip also contains two 
cross structures without antidots ("plain" sample). 

The fivefold Penrose lattice consists of two types of 
rhombuses with equal sides ap = 1.54 /im: "large" and 
"small" ones, with angles (20, 36) and (6, 46), respec- 
tively (6 = 36°). Accordingly, the short diagonal of the 
small rhombuses ap/r = 0.952 /xm; r = (1 + v5)/2 « 
1.618 is the golden mean. The rhombuses have been ar- 
ranged using inflation rules in order to assure the 
quasiperiodicity of the lattice. In the triangular lattice 
the next-neighbor distance is a<r = 1.49 /im. For the 
random arrangement of antidots their (x, y)-coordinates 
were generated by a 2D array filled with uniform random 
numbers and then scaled to give n p = 0.52 /im -2 . 

To characterize our devices, we measured resistance R 
vs. T (at B — 0) and determined T c of the different 
bridges on each chip. The left inset in Fig. Ufa) shows 
R(T) curves for perforated bridges and for one bridge 
without antidots on chip #2. For the perforated bridges, 
T c = 8.660 K is reduced by 12 mK compared to the bridge 
without antidots. This can be attributed to a small con- 
tamination from the resist structure during Nb film depo- 
sition, as we observed a similar behavior on other chips. 
The perforated bridges also show a larger normal resis- 
tance, which can be ascribed to geometry effects and to 
a slight reduction in the mean free path £ as compared 
to the plain bridge; again, this is consistently observed 
for all chips. With this respect, the R{T) curves shown 
in the left inset of Fig. |2ta) are representative for all de- 
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FIG. 2: (Color online) Comparison of four different bridges 
(chip#2, r = 200 nm). Main graphs: I c vs B (V c = 2/jV); 
T decreases from (a) to (c). Insets in (a): R vs T at B — 0, 
/ = 10/xA (left) and I c vs (1 - T/T c ) at B = Si/2 (right); 
grey solid line is calculated with a simple core pinning model. 
Insets in (b),(c): / vs V at B = (right) and B = B\ (left). 



vices on various chips which we investigated. From the 
measured resistivity Piok = 5.52 /xOcm of the plain film 
at T = 10 K and the relation pi = 3.72 • 10~ 6 /illcm 2 H2 
we estimate £ « 6 nm. In total we investigated 12 bridges 
on three different chips (#1, #2, #3). Below we present 
data for bridges with r = 200 nm antidots on chips #2 
and #3. Bridges with r = 125 nm antidots on chip #1 
behaved similarly. 

Figure |31 shows I c (£?)-patterns of a Penrose array at 
three different temperatures close to T c . Here, I c is de- 
fined by a dynamic criterion, i.e., via the detection of 
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B/B 1 

FIG. 3: (Color online) I c vs B of a Penrose array (chip#3, 
r=200 nm, V c = 1 ^V) at different T close to T c = 8.425 K. J c 
is normalized to I C (B = 0) = 0.71mA, 0.39mA and 0.16mA 
at T/T c =0.9976, 0.9982 and 0.9988, respectively. Inset shows 
magnification at 0.5SB/BiS1.15. B has been swept through 
a full cycle (from its minimum value up to its maximum 
and back). Dotted lines show simulation result redrawn from 
Fig. 2(e) in O for N p = 301. 



a threshold voltage V c — 1 /xV and the T-stability is 
S3 1 mK. For comparison, we insert the calculated I C (B) 
dependence (as dashed line), replotted from Fig. 2(e) in 
[3 (for A^=301 pinning sites), without any adjustable 
parameter. For these calculations a static I c definition 
(number of pinned vortices over the number of vortices) 
was used. However, Misko et al. 13] mention that 
dynamic simulations give essentially equal results. We 
find very nice agreement with our experimental data for 
the two highest temperatures (lower and middle traces): 
The narrow peak in I C (B) at B S3 broadens above 
B/B\Ki 0.2, until a local minimum in I C (B) is reached at 
B/Bi S3 0.5; a two-peak structure appears with I c max- 
ima at B/Bi — 0.81 w r/2 (broader peak) and at B = B\ 
(narrow peak). Above B\, I c drops rapidly with increas- 
ing B. The two-peak structure was predicted in [13J for 
the case that the gain in pinning energy E p i n exceeds the 
increase in elastic energy E e i associated with the defor- 
mation of a triangular vortex lattice. Misko et al. p~3| 
find the broader peak at a value of B v / t /Bi = 0.757, 
which corresponds to filling of only three out of four 
of the pinning sites on the vertices of the small rhom- 
buses, except for the N4 cases when two small rhombuses 
are connected along one side, and share one such va- 
cancy. We find for our large Penrose arrays Nd = N s /t 2 
(N s and Ni are the number of small and large rhom- 
buses, respectively). With the number of vertices in 
the Penrose lattice N p = tNi — t 2 N s , one obtains 
By/t/Bi = 2/t 2 S3 0.764. This is just slightly below the 
matching peak which we observe at r/2 = (N p — ^±)/N p . 
The latter corresponds, e.g., to a vortex configuration 
with every second small rhombus having only three out 
of four pinning sites occupied. In addition, we find at 
the highest temperature T/T c — 0.9988 a third peak in 



I C (B) at B/Bi = 0.65 ~ 1/r, which can be associated 
with filling of only three out of four pinning sites on the 
vertices of all small rhombuses. 

With decreasing T another peak-like structure appears 
at B/Bi = 1.77, close to 1 + (t/2). In the same T-range, 
the triangular array with same antidot size (on the same 
chip) develops a clear peak in I C (B) at the second match- 
ing field Bi- This can be explained by an increase of the 
saturation number [23l| n s S3 r/2£(T) above a value of two 
(£ is the coherence length). I.e., by lowering T, pinning 
of multiquanta occurs. The observed matching peak at 
1 + (t/2) corresponds then, e.g., to double occupancy of 
antidots, except for one antidot on every second small 
rhombus, which is only occupied by a single vortex. The 
fact that this matching peak is not very pronounced in- 
dicates that E P i n only slightly exceeds E e i, which also 
explains the missing of a matching peak at B = B 2 ■ 

The matching peaks in Fig. [3] are most pronounced 
at the highest temperature, T/T c = 0.9988, and gradu- 
ally transform into a plateau-like I C {B) pattern (within 
0.5 < B/Bi < f) at T/T c = 0.9976. This plateau is a 
remarkable feature, which we associate with the effective- 
ness of the many built-in periods of the Penrose lattice. 
We note that the Pearl length [H A(T) tx (1 - T/Tc)- 1 , 
which sets the vortex interaction range, changes from 
S3 80 to 40 /xm with decreasing T in Fig.0 In contrast to 
the calculations in [Tflj , with a penetration depth of the 
order of the pinning lattice spacing a, we find distinct 
matching peaks only for A»a, i.e. if the vortices inter- 
act over a large fraction of the Penrose array. Certainly, 
a better understanding of our results requires a more de- 
tailed theoretical analysis, including the T-dependence of 
the vortex interaction range, pinning strength and pin- 
ning range, as well as thermal fluctuations, which might 
be important, in particular close to T c . 

For direct comparison, we simultaneously measured on 
a single chip (#2) three perforated bridges with differ- 
ent antidot configurations (Penrose, triangular and ran- 
dom; r = 200 nm) and one plain bridge. Figure [21 shows 
I C (B) at three different values of T/T c . In contrast to the 
'plain' and 'random' samples the Penrose and triangular 
arrays show clear matching effects, with identical B\, as 
designed. The triangular array shows very pronounced 
matching peaks at B\ and B2 . We do not observe higher- 
order matching peaks, which indicates that n s S3 2. 

The I C (B) pattern for the Penrose array is very similar 
to the one on chip#3 with same antidot size (c.f. Fig.|2J). 
Decreasing T, the multiple peak structure at B < B\ 
turns into a broad shoulder. When T is lowered further 
[see Fig. Efc)], the shoulder transforms into a dome-like 
structure. We cannot give a concise explanation for the 
shape of this very broad central peak; however, we note 
that we observed this on all investigated Penrose arrays. 

Comparing absolute values of I c for the Penrose and 
triangular array shows that very close to T c critical cur- 
rents at B = 0, B\ and Bi are quite similar; however, 
due to the stronger reduction in I c of the triangular array 
between matching fields, the Penrose array is superior, in 
particular for small fields, below B\ [see Fig.^a)]. This 
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situation changes with decreasing T, as the reduction in 
I c of the triangular array between matching fields be- 
comes much less pronounced [see Fig. ^c)]. I. e., we 
cannot confirm the prediction |l3j | that the Penrose lat- 
tice provides an enormous enhancement in I C (B) over the 
triangular one, except for T very close to T c and fields 
between B = and the first matching field. This is 
also visible in the right inset in Fig. Ufa), which shows 
I c (l — tjt) at B = 0.5£?i for all four samples. Above 

1 — 3 • 10~ 3 the Penrose and triangular array 

show almost the same I C (T), which approaches, within 
a factor of two, the Ginzburg-Landau depairing current 
Id p oc iJ 2 £ oc (1 — y-) 3 / 2 , with the thermodynamic crit- 
ical field H c . This scaling can be derived within a sim- 
ple core pinning model, neglecting vortex-vortex inter- 
actions. The resulting I C (T) - shown as grey line in the 
right inset in Fig.|2ta) - is then simply determined by the 
maximum pinning force of an antidot, which we calcu- 
lated, following the approach in Q based on the London 
approximation |23j |. The 'random' sample gives always 
significantly smaller I c as the triangular and Penrose 
arrays. However, it shows a steeper / c (T)-dependence, 
which scales as (1 — ^r) q over the investigated range of 
T, with an increase in q from g = 2ati? = Oto<7 = 2.2 
at B = O.bBi to q = 2.4 at B = B\. This behavior can- 
not be explained within our simple model and certainly 
deserves further investigations. 

Finally, we note that the choice of the voltage crite- 
rion V c for determining I c does not significantly affect 
the shape of the / c (B)-curves; we checked this for various 
samples and V c ranging from fOnV to lOfiV. For a more 
detailed analysis, we measured current voltage character- 
istics (IVCs) as shown in the insets of Fig. |2Ib),(c) for 
all four bridges at T/T c =0.9977 and 0.9954, respectively. 
The IVCs of the perforated bridges evolve smoothly into 
the resistive state. In particular, there is no sudden 
change in differential resistance Rd pointing to a sud- 
den increase of the number of moving vortices. At B = 
(right hand insets) the IVCs of the different bridges do 



not cross (this holds up to at least 1 mV). Such a cross- 
ing is, however, observed in finite fields. As shown in 
the left hand insets (B — B{), Rd of the plain bridge is 
smaller than for the other bridges. At V « 300 /iV its 
IVC intersects those of the triangular and Penrose ar- 
rays. A similar crossing is also observed for other values 
of field and temperatures, however, always well above the 
threshold voltage V c used for determining I c . 

In conclusion, we experimentally verified theoretical 
predictions of the main features of the critical current 
dependence on the applied magnetic field in a supercon- 
ducting thin film with a quasiperiodic Penrose lattice of 
antidots at temperatures close to T c . In particular, we 
found matching of the vortex lattice with the quasiperi- 
odic lattice of pinning sites very close to T c , and we asso- 
ciate various matching peaks in I C (B) with distinct reg- 
ular arrangements of the vortices. In addition, we di- 
rectly compared different arrangements of artificial pin- 
ning sites in our Nb films. We find an enhancement of I c 
in films with Penrose lattices as compared to films with 
random arrangement of pinning sites and films without 
artificial pinning sites, however no significant enhance- 
ment of I c as compared to films with triangular antidot 
lattices. With respect to applications, it will be interest- 
ing to perform more detailed investigations on the effect 
of optimum antidot size (i.e. saturation number) and den- 
sity over a wide range of temperatures in quasiperiodic 
pinning arrays. 
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Note added. After first submission of our manuscript 
we learned about related work by Villegas et a/.[2{| and 
Silhanek et aL[2(| on quasiperiodic pinning arrays in su- 
perconducting films with magnetic dots. 
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